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As part of the Ministry of the Environment's hydro^ 
logic research studies , the River Basin Research Section is 
developing mathematical models to describe the complex 
interactions of the various components of the hydrologic 
cycle in research watersheds. During model development, it 
is useful to apply developed techniques to existing hydro- 
logic questions in order to test the viability of these 
techniques for solving problems related to local hydrologic 
conditions . 

The results of the application of the streamflow genera- 
tor described in this study will form part of an overall 
comprehensive report on the Thames River basin management 
study, to be released at a later date. The development and 
application of this daily streamflow simulator is outlined 
separately here to indicate the use of such a technique 
possibly as an aid to others engaged in the simulation of 
streamflows . 




K. E. Symons , Director, 

Water Quantity Management Branch 



Toronto, January 1, 19 74. 
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FOREWORD 



In generating stochastic sequences of flow for the 
Thames River stations, using statistical properties of 
historic data, it is extremely important to realize that 
the sequences are generated by statistical methods that 
do not pretend to provide cause-effect models for actual 
flows. Fiering and Jackson (3) indicate that the data 
generated by using such methods "while neither actual 
historic records nor predictions of future flows, are 
close enough to possible (but not observed) historical 
records that may be used to determine, within defined 
statistical errors, the several quantities of interest," 
from the point of view of a water management study. 
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DAILY STREAMFLOW SIMULATION FOR 



THE THAMES RIVER BASIN 



INTRODUCTION 



In participating in the Thames River basin water 
management study, the Water Quantity Management Branch is 
presently undertaking several studies dealing with the 
collection, compilation and. analysis of gro'und and surface 
water data. In order to provide essential input to the 
water quality simulation model under development by the 
VJater Quality Branch, the Water Quantity Management 
Branch was requested to provide extended daily strecmiflow 
data sequences for several stations on the Thames River. 

The objective of the study was to derive a suitable 

multiple-station daily streamflow generator for the 
Adelaide, I'aling and Byron stations in the vicinity of the 
City of London, using an acceptable mean monthly generator 
to provide the proper required input of synthetic monthly 
means to the daily flow generator. 
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THE THAMES RIVER BASIN 

The Thames River, the headwaters of which are located 
in the highlands of the counties of Perth and Oxford, 
northeast of London, drains approximately 2,250 square 
miles. The river carries water for more than 190 miles 
to Lake St. Clair. The Thames River above London (Upper 
Thames) has numerous tributaries, all of which drain into 
two main branches: the North Branch (app. 660 sq. miles) 
and the South Branch (app. 5 30 sq. miles) , the confluence 
being within the City of London. Below London., the Thames 
River (Lower Thames) flows through a comparatively narrow 
drainage area where all the tributaries appear relatively- 
small (9) (see Figure 1, Appendix II). 

'^lic mean annual temperature in the basin ranges from 
4 8.7°F at Chatham, near Lake St, Clair, to 43.o^F at 
Mitchell in the Upper Thames. The mean annual precipitation 
increases generally in an upstream direction from 31.59 
inches at Chatham to 39.01 inches at Stratford, 

fites Used in the Study 

The daily streamflow generation is to be applied 
simultaneously to three sites (Figure 1) : 

1. Adelaide station (02GE002) on the North Branch 
Thames River 

2. ialing station (02GD001) on the South Branch Thames River 

3. Byron station (0 2GD00 3) below the confluence of the North 
and South branches (Thames River) . 

The historical records available for the three stations are 
in the form of sequences of discrete values of average daily 
streamflows at each gauging station. Table 1 in Appendix I 
gives a summary of available data for these stations. 

The eight-vear period 1923-19 30 was chosen as the 
concurrent period of record for the three stations , the 
historical data for which represent unregulated flow sequences « 
It was recognized that a sample of streamflow data of this 
size is a small one; however, a decision was made to attempt 
to generate synthetic data on the basis of the above-stated 
availability of data and for the additional two reasons: 

i - The available concurrent records of regulated flows 
cannot be used because the extent of regulation has 
not been well-docum.ented and therefore the effects 
cannot be defined. 

ii - 1^.1 >i comparisons of the means and standard deviations of 
annual flows for the eight-^^ear period to those for 
longer, non-concurrent, unregulated periods available 
for stations 02GD001 and 02GD00 3 (Table 2) were within 
an acceptable range and no adjustments were considered 
necessary to the parameters for the shorter period. 
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LITERATURE REVIEW OF DAILY STREAMFLOW SIMULATION TECHNIQUES 
General 

In hydrologic studies, river discharge, as a variable 
observed or recorded in time, can be considered as an 
unidimensional (time), ^.tochastic hydrologic process. 
Sequences of streamflow data characterized by statistical 
properties are called stochastic sequences (processeis) and 
as they are observed or recorded in time, they constitute 
a time series (12) . Hydrologists are often handicapped 
because of the short sequences of available flow data. As 
a solution to this problem, hydrologists during the past 
decade have devised various data simulation techniques 
which may be applied to historical records to obtain 
hypothetical sequences of data which preserve the significant 
statistical properties of the historical data. All these 
simulation techniques are mathematical models that make 
use of the computer and programming logic. In this sense, 
simulation is a tool that combines theory, data, and programming 
logic to express in mathematical terms the pertinent elements 
of a complex hydrologic problem (10) . 

Daily Streamflow Simulation Models 

In the past, stochastic models have been applied to 
annual and monthly stream.flow data, obtaining results 
compatible with historical flows (Thomas and Fiering (8) , 
Beard (1), and others). Extensions of these models for 
generating daily flows were first attempted by Halter and 
Miller (5) , who developed a daily flow" simulator using a 
linear regression model which generated 30 flows each 
month, using as its basis the mean monthly flow and the 
standard error of the monthly flow. The simulated hydro- 
graphs, nowever, \7ere not adequate because the serial 
correlation among previous flows was not incorporated into 
the generator; however, flows during recession periods were 
considered to be exponentially correlated. 

At the International Hydrology Symposium (Fort Collins, 
1967), Beard (2) presented a single-itation model based on 
generation of monthly streamflows and subsequent allocation 
of the monthly total amounts to each day. Ilie monthly 
streamflow generator used by Beard was that developed in 
the Hydrologic Engineering Centre of the Corps of Engineers 
(4). Beard's miodel thereby provided a realistic representation 
of daily flows which met the four following specifications: 



-. ^ _ 



i - the generated flows must maintain the expected mean 
of the observed flows , 

ii - the variance about the monthly mean must be comparable 
to that of the historical record, 

iii- discontinuity between daily flows at the end of one 
month and the beginning of the next month must be 
avoided, 

iv - the generated flows must be properly related in time (7^ . 

Payne , Neuman and Kerri (6) , subsequently introduced 
a daily streamflow generation suitable for an arid region, 
in which they suggested the re -arrangement of the times of 
occurrences of flood peaks within an annual cycle, in 
order to preserve the pertinent properties of the basin 
hydrograph. In all other respects, the Payne et al model 
was the same as Beard's model. The operational device as 
suggested by Payne et al, however, appears unsuitable for 
humid areas where there are usually two to five peaks per 
month during the rainy season. Under such varying conditions , 
the re-arramgement of the historical record would likely 
pose a serious problem. The re -arrangement of the flow 
peaks to occur at one or two specific dates within one 
month, or a season, destroys the time sequence, and thereby 
reduces the accuracy of the response of an existing system 
to potential non-historical flow sequences, because the 
intervals between consecutive peaks or lows has been 
artificially fixed. 

In 19 72, Environment Can ada (7), introduced a multi- 
site daily streamflow generator for the Fraser River basin. 
Tliio generator used basically the same technique as described 
by Beard and incorporated the following two additional 
specifications for a multi-site model: 

i - pynthetic daily flows at different sites are properly 
related to one another in terms of location, 

ii - t ' e hydrologic lag-time between upstream and downstream 
stations is preserved. 

The Fraser River model also uses synthetic monthly streamflow 
as an input. These are generated using the technique as 
developed by Young and Pisano (11) . 

It should be noted that, based on observed, historical 
records, statistical relationships can be derived for daily 
flows to allow the re-generation of the observed record; 
however, in order to generate sequences of daily streamflows 
for periods greater than the historical record, synthetic 
monthly means covering the same period for which daily flows 
are to be generated are required, as input to the daily 
generation program. 



In reviewing the various techniques of simulation, 
it was concluded that Beard's model likely represents the 
current state-of-the-art. r.eard'f^. original model, however, 
was designed for a single-station problem. Iheref ore , 
beard's technique, as modified by Environment Canada for 
a multi-site study in the Fraser River basin, extended to 
suit the Ihames River, was considered to be most effective 
in achieving the objectives of the proposed study. 
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BEARD'S DAILY STREAMFLOW MODEL 



The daily streamflow generation as used by Beard 
consists of a second order Markov chain, using standardized 
variates. For Markovian models, such as Beard's, each 
value of the process at time t depends on its immediate, 
previous values at times t - 1, t - 2..,,.., and on an 
independent random component. In mathematical notation 
this can be written as follows: 



^t = 



^t-i + 



X 



t-2 



+ 



a X^ + E. 
n t-n t 



... (1) 



where 



X. = the present value of the process, 

X , ...X - the past values of the process. 



a,.,. a = coefficients, 

1 n ' 

E «= a random component with a certain 

probability distribution or pattern. 

The parameter n gives the order of the process (i.e. 
order 2 in Beard's model where the present flow value is 
affected by flow values for the two previous days) and 
also represents the measure of its memory. Memory applies 
to the span of time during which an event has meaningfully 
affected the events following. 

The generation of daily flows in Beard's model is 
accomplished by two computer programs, an analysis program 
and a generation program. 

In the analysis program, daily flows are first 
classified by calendar month, and their frequency char- 
acteristics (logarithmic mean, standard deviation, skew 
coefficient, serial correlation coefficient, and the 
correlation coefficient between flow logarithms of each day 
and the second antecedent day) for each calendar month are 
determined from observed data. 



Standardized variates are obtained by subtracting the 
mean logarithm from the daily flow logarithm, dividing 
this value by the standard deviation, and transforming to 
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normal by use of the following equation, /hich approximates 
the Pearson Type III distribution: 

t*= jg ({g i: + 1)^^^ - 1) + g ... (2) 

g ? 6 

where 

t*= a normal standard deviate, 

K = a Pearson Type III standard deviate, 

fj = the skew coefficient. 

The analysis program proceeds with computing the 
linear regression coefficients of the standard deviation 
of daily flow logarithms within each month of the record 
to logarithms of total flow for each month. 

The generation program of Beard '■,> model generat€;s for 
each day standardized variates conforming to the serial 
correlation observed in the recorded data for that calendar 
month, using the following equation: 

where 

X. ,X. ^ ^^'+2 ~ standardized variates for successive days, 

X = a random standardized variate, 

Y 
b, , b^ - regression coefficients, 

=2 
R = the determination coefficient for the 

regression equation. 

The generation program obtains the logarithm of the 
monthly mean flow as a first estimate of the mean logarithm 
of daily flows to be generated. It is noted that the 
generation program requires a set of monthly means 
(logarithmically transformed) to initiate the generation of 
daily sequences. J-^ the available monthly means of the 
historical record are used as input to the generation program, 
along with the frequency characteristics of the same record 
obtained from the analysis program, then the generation 
program will reproduce (theoretically) the historical daily 
sequence. This characteristic of the program is used to 
check the ability of the model to reproduce the historic 
hydrograph of a particular year or years of the record for 
visual comparison. 






In order to generate sequences of daily stream flows 
that are longer than the available observed record, but 
are statistically similar to the historical record, a 
set of synthetic monthly means of the same length is 
needed as an input to the generation program. As was 
stated earlier. Beard '3 model utilizes the monthly stream- 
flow generator developed in the Hydrologic Engineering 
Centre of the Corps of Engineers, to obtain such synthetic 
monthly flow sequences, 

'The daily generation program transforms, the standardized 
variates to logarithms of daily flows by use of the following 
approximate transform equation for the Pearson Type III 
distribution: 

K ^ 2 (g (t* - g) + 1)-^ - 2 ... (4) 

g t t g 

which is simply another form of Equation (2). The logarithms 
are converted to flows which are added to obtain a monthly 
total. The sum includes random components and is therefore 
somewhat different from the desired total entered from the 
monthly streamflow model. A ratio of these totals could be 
applied to each daily streamflow, but this would ordinarily 
result in discontinuities at the end of each month. In 
order to smooth these discontinuities, Beard' '^. model preserves 
the random numbers used for that month and generates new 
daily values o^n 'th^e basis of a temporary monthly total. 
The latter is obtained by multiplying the desired monthly 
total by the ratio of the desired total to the generated 
total, for the first pass. 

The flows resulting from the second pass generation 
are very nearly the total of the desired monthly amount. 
These flows are therefore multiplied by a ratio to obtain 

the exact monthly total, so that a serious discontinuity 
at the end of each month is not created. 
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GLuLl^Vi'ION OF MONTHLY STREAMFLOWS FOR THE THAMES STATIONS 
BASED ON THE YOUNG AND PISANO MODEL 

The daily streamflow generator, as originally used 
for the Eraser River study, is dependent on the existence 
of an acceptable monthly mean generation model to provide 
the proper input. Tie mathematical technique, used for the 
generating of the necessary monthly streamflows for the 
Thames stations , is based on the Young and Pisano meithod 
(11) , Ihe technique appears more suitable than others 
in dealing with short periods of record for generating 
synthetic sequences , as all months are combined for the 
correlation analysis, rather than treated independently (7). 

The method presumes an underlying Markovian structure 
and considers cross -linkages or interrelationships between 
stations along both the time and space axes. The advantages 
of the method are its straightforward approach to the 
problem of flow generation, its computational ease and wide 
applicability . 

The objective, in general terms, is to extract from the 
historical records , statistical parameter values to construct 
a model of the monthly historical time sequence, in the 
following manner: 

Consider Y years of data at n sites. The data constitute 
an observation matrix X, of monthly values, having dimension 

n . t (t = 12Y) . The matrix X is used to construct the 
model which generates S years of synthetic data. 

rollowing the Young and Pisano principles , the computer 
program for this model was written by the Engineering Division, 
Inland Waters Directorate, Department of the Flnvironment . 
The program proceeds along the following steps : 

1, Computes the residuals,, using the following equations 

V 

r. /' X . . K- m^^ , /r- \ 

13 = 13 Kj .... ( 5 ) 

"Kj 



where 



r. . - are the residuals, 

X. .^ "' are individual elements of matrix X, 

13 

S^. ^ are the corresponding standard deviations, 

m , "* is the average monthly value (an average of the 



I-D 



monthly averages). 
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K - denotes the s ite {K - 1 , 2 , , . , ,n) , 
i - denotes the year (i = 1, 2,,..,Y), 
j - denotes the month (j - 1, 2,... ^12), 

Each element in matrix X is replaced by its corresponding 
residual to yield a residual matrix R. of dimension 
n . t (t = 12 Y) , 

Transforms the historical record in matrix X, by taking 



login X. .K /x, . of each element and transforms 

^10 1] ' ij 

further the new matrices into R^ , R^ matrices. 

3. Compares the matrices R^ , R„ , R- and adopts the R matrix 

with the smallest average skewness. As normal data have 
an expected skewness of zero, the strategy is to pick up 
the option which yields a minimum average skewness. 
Average skewness is the sum of skew coefficients at 
each site (row) divided by n. Ihis initial operation 
was referred to by Young and Pisano as the finding of 
the MST (Minimum Skewness Transformation), 

4. Computes the means and standard deviations of the 
normalized data, r^atrix R, for each calendar month of 
each site and then standardizes the data. At this 
stage the data are assumed to be normally distributed 
with zero mean and un,it standard deviation . 

5. Computes the matrices Mo and M_, which are the sample 

estimates of the variance-covariance matrix and the 
lag-one covariance matrix of the residuals, using the 
following equations : 

T 
' '^' ..,(6) 

... (7) 



Mo = RR 


(1) 




n 


— T 
M_ = RR- 


(1) 
n 



where 



Mo - is the variance-covariance matrix (the cross - 
correlation matrix) , 

M_ - is the lag -one covariance matrix (the lag -one 
1 correlation matrix) , 

R - is the residual matrix, 

R - is the same as R except that the rows have been 
shifted one time frame and the_ first column of 

R becomes the last column of R, 

m rp 

R , R^ are transposed R and R matrices. 
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Sets up the generation model which is based on the 
following equation: 



where 



R„ ■- is the £ column of residuals, 

Ai;, -• is an n x n matrix, 

B - is an n X n matrix, 

e — is a n X 1 vector of independent (0,1) normal deviates, 

g - is the time frame . 

The matrix A is obtained using the following equation: 

A = M T M "-^ • .. (9) 

-1 o * 



where 



M_ — is the lag-one covariance matrix, 

-1 . . 

M — is the inverse of M matrix. 
o o 

The m.atrix B is obtained using the following equation! 

BB*^ = M - M • M ""^ • m"^ T ... (10) 

O ~1 o -X 

where 

T 
1 - is the transposed B matrix, 

I'i"^__-] - is the transposed M , matrix, and all otlier terms 
are as previously de fined, 

llie remaining analysis pertains to the solution of 
Equation (10) . The matrix B is obtained by assuming 
B to be a lower triangular matrix. This assumption 
allows a solution for the coefficients of B. 

fm Generates synthetic residuals: (1) set R = (all 

o 

elements); (2) calculates R, and succeeding values, R^ , 
from Equation (8N 
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3. Destandardizes the data using the following equation: 
X ^ = S,,. • r. .^ + m,, . . ,. (11) 

where all terms are as defined for Equation (5), 

9. uses the inverse of the adopted normality transform to 
produce sequences of synthetic monthly data. 
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GEKERATION OF DAILY STREMIFLOWS FOR THE THAMES PIVER 
STATIONS BASED ON THF ERASER RIVER MODEL 

The mathematical techniques employed in the generation 
of daily streamflows for the Thames River stations are 
basically the same as those employed in the Fraser River 
basin study by the Engineering Division, rnvironment Canada. 

j le generation of daily streamflows is accomplished by 
two computer programs. The first program, is restricted to 
an analysis of the historical daily streamflows and the 
computation of statistical parameters of the histo^rical 
data, an are required in the generation procedure. Th,e 
second computer program em.ploys these parameters to generate 
synthetic daily streamflow data for a specified number of 
years for which synthetic (from Young and Pisano monthly 
generation) or historic monthly flows are available. The 
synthetic daily flows were generated for the calendar 
months of May to November (inclusive) , at three sites 
simultaneously (the Adelaide, Ealing and Byron stations). 

The Historical Data Analysis Program 

The program reads the station identification and 
location information for the three stations in question: 
(Byron-02GE002, Ealing-12GD001 and Adelaide-02GD003) ; it 
also reads the eight years of concurrent historical data 
(192 3-19 30) for each. The data analysis is performed for 
one month at a time and proceeds as follows : 

1. Computes the monthly means (K) , standard deviations (S) 

and the monthly totals for the monthly flows , 

n 

X - i- . f , X. , (i - 1, 2,3,,. .n) . .. (12) 

n X — 1. X 



n 

E 

'n - 1 i = 1 "i 



S — i— ^ . I ^ X.^ - ^_n_ X^, (i = 1, 2. . .n) ... (13) 

n-1 



2. Transforms daily flows to log^-. values. 

3. Computes the monthly means, totals and standard 
deviations of transformed data, 

4. Regresses the standard deviations of the transformed 

data O'n their monthly tO'tals for each station and 
stores the regression coefficients and constants.. 
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computes long-t/irn means , standard deviations and 
skew coefficients (y) by month for the transformed 
data f 



n n 

13—12 —3 



y ^ 



n 1 I X. 
n 1=1 



^ 2 .-2,1.5 



(i = 1, 2, 3. . .n) ... (14) 



6. Standardizes the transformed data by subtracting the 
mean logarithms from the daily flow logarithms and 
dividing the result by the standard deviation. (The 
standardized variates were considered to be log- 
normally distributed and the normalization step used 
in Beard's model was omitted from the Thames River 
application) . 

7. Performs a regression analysis for each gauging station. 
Daily data are regressed on flow in two antecedent days 
and on the lagged flow at the doivnstream station, 02 GEO 02. 

8. punches out statistical parameters of the historical 
data to be transmitted to the data generation program. 

Input Data 

Card input data to the analysis program are in the 
following order: 

1. Gauging station information cards. 

2. Magnetic tape input consisting of daily streamflow data 
at the three sites: Byron, Adelaide and Ealing. Data 
are in the Water Survey of Canada format 6 7-002. 

Output Data 

The historical data analysis program prints out and 
punches the following information: 

1» Standard deviations and long-term means of May mean 

flows in log, Q form for the three gauging stations, (initializing 
data) . 

r . Gauging station information. 
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3. Rogression constants for the three gauging stations 

derived from the regression analysis between historical 
monthly total flows and the standard deviations of 
daily data, for each month, about the monthly mean. 

4* joegression coefficients for the three stations derived 
in the regression analysis outlined under item 3, 

3, Regression coefficients for the seven months and three 
stations derived in the regression analysis between 
daily flows at the dependent (Byron) and independent 
stations (Ealing and Adelaide) . 

6. Kegression constants for the seven months for the three 
stations derived in the regression analysis outlined 
under item 5 . 

7. The standard errors derived in the regression relation- 
ship between daily flows at the dependent and the 
independent stations • 

8. Average standard deviations in log, « units of daily 

flows about the monthly mean for each month at ectch 
station, is derived from the historical record. 

Daily Streamflow Data Generation Program 

The generation program produces a specified number of 
years of streamflow hydrographs over the required period of 
May to November, at the three gauging stations in the Thames 
River basin. Ifiese flow sequences are based on the 
statistical information derived from the historical daily 
streamflows at each of the three stations. The generation 
program accepts the statistical parameters prepared during 
the analysis of historical data, and a set of synthesized 
monthly means (see generation of synthetic monthly means) 
and proceeds as follows : 

1* Reads the average standard deviations of daily flows 

about the monthly means for each station in log, ^ units. 

2. Reads the monthly mean values as a first estimate of the 
mean daily flows to be generated for all stations being 
considered . 

$, Computes the monthly totals using the following equation 

monthly total = mean X number of days in month. 



4. Transforms the monthly totals and means to log^^ values. 

5# Defines the flow values for the two antecedent days 
(2 9th and 30th of April, in this case) to initiate 
annual flow sequences . The starting values for May 
are set as May monthly means for that year, standardized 
and normalized with respect to the long-term statistics 
for the month, 

6. Defines the total lag for the downstream station used 
in the regression analysis, 

7, Sets the monthly standard deviation at each location 

to be equal to the average value of the monthly standard 
deviations obtained in the analysis program. 

C. Generates for each day of each month, flow sequences in 
the form of standardized variates using the following 
equation: 

^i+2,Yr ^1 ^i+l,n ^ ^2 ^i,n ^ ^3 ^i+2 + il-^,n-l ^ 

^4 ^i+2+£^,n-2 + ^1-R^ .X ...(15) 

<i r 

where 

X. , X. . T , -^-.o ~ cire standardized variates for 
i.n 1+1 ,n 1+2, n . j , . .• 

' ' successive days at station n, 

X X 

i+2+il-,n-l, i+2+5,2 ^"2 ~ are standardized variates at the 

' stations downstream of station n 

on day i + 2 + £, where £], represents 

the lag time between n and n-1, 

and ilo the lag time between n and 

n-2 (the lag time for the Thames 
River study was assumed to be equal 
to zero because of the proximity of the 
_g, stations to each other) , 

m - is the determination coefficient for the regression 

equation. 



b 



^, b^, b_, b. - are regression coefficients, 



X - i.s a random standardized variate. 
r 



'^ 
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The sinilarity between Equation {15) and Equation 
(3) in Beard'-J r^odel is obvious, the only difference 
being that Equation (15) takes into consideration 

the effect of the flow at stations n-1, n-2 on the 
flow at station n. 

Transforms the sequences of standardized variates 
into logarithms of daily flows by multiplying each 
standardized variate by the log value of the standard 
deviation and adding to that the log value of the mean* 
lliese logarithms are then converted to naturalized flows, 
which are added to obtain a monthly total. 

'.the final monthly total is computed in the program 
through two passes (iterations) as described in Beard's 
model, 'hen the first iteration is completed, a new 
target total is set for the second iteration. After 
the second iteration is completed, the program performs 
a correction by multiplying each generated daily 
streamflow by a correction factor. The correction factor 
equals the ratio between the desired monthly mean and 
the generated monthly mean inputs • 

10. Checks on month-end discontinuities. Tf the discontinuity 
is v/ithin an acceptable limit, the sequence of generated 
data for a station and a month is accepted and stored. 
Cthen^ise, (i.e., if discontinuity is outside the 

acceptable limits) the flow sequences are regenerated 
up to ten times and the sequence with zero or minimum 
discontinuity is accepted and stored. 

Input Data 

The input data for the generation program are the 
statistical parameters derived in the historical data 
analysis program and a set of monthly means for the irequired 
number of years to be synthesized, n- outlined previously. 

Output Data 

The program output consists of synthetic daily stream- 
flows which are printed O'Ut and written in card-image form,, 
three card images per month, on magnetic tape in the following 
format: 



_ 1 p _ 



Column Format Rema rk s 

Card 11 Al a letter "S" signifying synthetic 

data 
3-3 I1A2I3 a Water Survey of Canada 

identification code 
9-11 T3 year identification 
12-13 12 month identification (5 to 11 

for May to September) 
14 II the integer "1" denoting the 

first data card for the month 
15-75 1016 synthetic flows for the first 

ten days in the month 
79-80 16 the total number of daily 

flow observations in the 
month 

Card 2 



1 


Al 


as 


on card 


1 


3-8 


I1A2I3 


as 


on card 


1 


9-11 


T3 


as 


on card 


1 


12-13 


12 


as 


on card 


1 


14 


Jiill! 


an 


integer 


If "■ II 



c^.enoting the 
second data card for the month 
15-7j 1016 synthetic flows for days 11 

to 20 in the month 

Card 3 1 Al as on card 1 and 2 

as on card 1 and 2 



1 


Al 


3-8 


I1A2I3 


9-11 


13 


12-13 


12 


14 


11 



as on card 1 and 2 
as on card 1 and 2 
the integer "3" denoting 
third data card for the month 
15-CO 1116 synthetic flows for days 21 

to the end of the month 
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TESTING OF THE GENERATED DATA 



General Remarks 



In synthetic data production, the generated data must 
preserve the significant statistical properties of the 
historical data. As Fiering and Jackson (3) indicate, the 
most that any generating scheme can promise about any 
statistic (mean, standard deviation, etc.) of the historical 
data is that the expected value of the same statistic is 
the desired one. The expected value of a statistic is 
its average value in an infinitely long sequence. However, 
because the generation algorithms are only used to form 
finite sequences of flows, the sample statistics obtained 
cannot be expected to be exactly equal to the historical 
statistics. 'Tliey O'nly tend tO' be near the historical 
statistics, and the closeness is expected to improve with 
the length of the generated data. Therefore, the preservation 
of the significant statistical properties of the historical 
record within the synthetic sequence implies that th,e 
expected values of these statistics are the specified values . 
The actual obtained value of a statistic for a finite 
synthetic sequence will be within theoretically anticipated 
sampling errors of the expected value. Hiis qualification 
should be taken into consideration during the examineition 
of the test results that provide comparisons between the 
means , standard deviations and coefficients of skew for 
individual months of the historical and synthetic data. 

Results of the Synthetic Monthly Flow Data 

h sequence of 1,00 years of monthly means was generated 
and O'ne set of 50 years, exhibiting the most extreme values 
in the low flow range, was chosen for the study. Low flow 
extremes were considered to be of particular importance in 
evaluating possible critical water quality conditions in 
the river. 

A series of tests were undertaken to check the resem- 
blance between the monthly historic record and that of the 
selected 50 years of synthetic data. These tests included 
the following: 

1. comparison between the means, standard deviations and 
coefficients of skewness of the historic record and 
synthetic data, 

2. comparison between the cross -correlation matrices and 
the lag-one correlation matrices, 

3. comparison between the auto-correlation coefficients, 

4. comparison between the cyclical characteristics of both 
the historic and synthetic data {spectral analysis) , 

5. comparison between "Hurst's K" coefficients. 
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6. comparison between the number of changes in quartile 
and the number of changes in the direction of movement 
for both the historic and synthetic data. 

Tables 3, 4 and 5 provide comparisons between the monthly 
means , the standard deviations of the monthly means about 
the long-term mean, and the average monthly skew coefficients 
of the historical flows and the selected set of synthetic 
flows. Corresponding values for both types of data in these 
tables are close but not identical. This is to be expected 
as the generation algorithm is designed to reproduce the 
historical statistics (in this case mainly the mean and 
standard deviation) only if an infinitely long sequence of 
synthetic data is generated. Different finite synthetic sets 
will have different statistics that are all distributed 
around the historical statistics within theoretically 
anticipated sampling errors. "^his prooerty allows the user 
to choose such synthetic sets that exhibit extreme values 
in evaluating, for example, the effects of floods or droughts 
for various water management studies. 

The comparison between the average means, standard 
deviations and coefficients of skew, as shown in Table 6 
for all months of the historic and synthetic data, show 
closer agreement indicating that the statistics for the 
finite synthetic sample tend to be near the expected values 
of the historic statistics. 

Table 7 shows the cross -correlation matrix Mo, the lag- 
one correlation matrix M_, , and the corresponding estimates 

of A and B matrices for the historical data, as indicated 
in equations 2, 3 and 4. The matrices Mo and M_^ of the 

synthetic flows are given in Table 8. Close agreement 
among corresponding entries in tables 7 and 8 is evident. 

Table 9 provides a comparison between the cross- 
correlation matrices of the historic and synthetic data 
on a monthly basis. 1"iese data also indicate general 
agreement. 

The auto-correlation test (Table 10) compares the 
correlograms of tlie historical and synthetic flows. As 
was stated earlier, the monthly generator presumes an 
underlying Markovian structure of the first order. This 
assumption implies that the flow during any month is 
considered by the model to be meaningfully affected by the 
flow in the previous month only. As it can be seen from 
Table 10, the correlograms of the historical record of lags 
higher than one have much smaller values than those for lag 
one. This shows that the presumed underlying Markovian 
structure of order one in the monthly generator is adequate. 
Table 10 indicates that the lag-one correlograms of the 
historical and synthetic flows are close. It also shows 
that the correlograms of higher order compare favourably. 

The cyclical characteristics of the historical and 
synthetic flows obtained by the spectral analysis test 
(Table 11) indicate general agreement. 
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rtroamflow data arc characterized by persistence. 
I.urst*!:. coefficient (K) is considered to be an indicator of 
the degree of persistence in streamflow data. Table 12 
compares the Hurst's K values for the historical and 
synthetic data. P.s can be seen from Table 12, K values for 
the historical and synthetic data are close, both being 
within the range 0.57 to 0.63, 

The tests for the number of changes in quartile and 
the number of changes in the direction of movement (tables 
13 and 14) for the historical and synthetic data indicate 
fairly good resemblance; liowever, the synthetic data show 
higher values. This is to be expected since 50 years of 
synthetic data are compared to only eight years of historical 
records . 

Results of the Synthetic Daily Flow Data 

ii 50-year sequence of syntJietic monthly means /generated 
accO'rding to the Yo^img and Pisano method^ v/'as used to generate 
a 50-vear sequence of daily flows. Two types of tests were 
undertaken to check the ability of the daily flow generator 
to reproduce realistic daily streamf lows : 

1. daily streamf lows were first generated based on available 
historical monthly means and the generated hydrographs 
were visually compared with the historical hydrographs 
for the period 192 3-19 30, 

2, a 50-year sequence of synthetic daily flows, as well 
as the historical record for the period 1923-1930 

(natural flows) and for the period 1956-1970 (regulated 
flows), '..;ere statistically tested and the results were 
compared on a montlily basis. In addition, the cumulative 
probability distributions of the historical and synthetic 
minimum daily flows for stations 02GD001 and 02GD003 were 
compared. 

Plots of the synthetic and historic hydrographs for 
the three stations under study are given in figures 2, 3 
and 4, for the year 1930, These plots illustrate the ability 
of the daily flow generator to reproduce realistic flow 
sequences . 

Table 15 provides a comparison between the means and 
average standard deviations of the historic record for the 
periods 1923-1930 and 1956-1970 and for 50 years of synthetic 
data on a monthly basis (May to November) , It is clear from 
Table 15 that the values of the above-mentioned two statistics 
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are very close for the historic period 1923-1930 (natural 
flows) and the period 1956-1970 (regulated flows) . This 
fact indicates that the regulation did not drastically 
change the characteristics of flow during the months of 
May to November. 'iable 15 also indicates that the means 
and average standard deviations for the synthetic data 
compare favourably with their corresponding values for the 
historical records. 

Table 16 provides a comparison between the skew 
coefficients of the above-mentioned three sets of data. 
Ihe comparison indicates that the synthetic data, in general, 
are slightly less skewed than the historical data. Table 
17, iowever, ohows that the long-term means and their 
standard deviations, as well as their skew coefficients, 
are in close agreement for both the synthetic and historical 
data. 

Table 18 provides a list of minimum daily discharges 
as recorded annually at stations 02GE002, 02GD001 and 02GD003 
for the period 1916-1970. In examining this table, it is 
evident that the extreme minimum daily discharges at 
stations 02GD001 and 02GD00 3 can reach as low as 2.0 cfs. 
'ihe regulated record at station 02GD003, for the period 
195 5-1970, contains flows as low as 0.5 cfs. Figures 5 
and 6 illustrate the cumulative probability distributions 
of historical and synthetic minimum daily flows for stations 
02GD001 and 02GD00 3. Extreme values for historical and 
synthetic flows are either close or identical in both figures, 
whereas the intermediate values are lower for the synthetic 
data. This indicates that this particular set of synthetic 
flows provides a lower envelope for probable flow sequences 
and consequently can give the user estimates of the expected 
frequency and indications of the severity of critical low 
flow conditions. 
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CONCLUSION 

A set of 1000 years of monthly streaimf lows , as well 
as one 50-year set of daily streamflows have been generated 
for the study. The statistical tests of the generated data 
indicate that the applied streamflow generator satisfactorily 
reproduced the most relevant characteristics of the historical 
records. 

The synthetic daily streamflows that were generated 
provide a lower envelope for probable flow sequences. 
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REC0I4MENDATI0N 

For a complete analysis of the sensitivity of various 
water quality control alternatives, 5t is advisable that 
additional 50-year sets of synthetic daily sequences be 
selected. This will allow the user to obtain all the 
possible flow ranges, including the medium and high ranges 
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TABLE 1 . Nummary of Available Data for Sites Used in the Study 



Station 
Name 



Station 
Number 



Drainage 
Area in 
mi^ 



Period of 
Record 



Remarks 



Thames R, 
at Byron 



02 GEO 02 



12 00 



Thames R, 
at Ealing 



Thames R. 
below 
Fanshawe 
(Adelaide) 



02GD001 



519 



02GD003 



560 



-continuous 
Oct, 1922- 
Sept.1931 

-spring flood 
discharges 
Sept. 1922- 
Aug.1955 

-continuous 
Aug, 19 55 - 
Dec. 19 70 

•continuous 
May 1915- 
Dec.1970 



■continuous 
June 1915- 
nept.1934 
Oct. 1935- 
Sept.1944 
Oct. 1955- 
Dec.1970 



records excellent. 
Since 1953 
discharges affected 
by regulation at 
Fanshawe Dam, nince 
1965 at Wildwood 
Dam, and since 196 7 
at Woodstock Dam 



records good. Since 
1967 discharges 
affected by regulation 
at Woodstock Dam 

records excellent. 
Since 195 3 discharges 
affected by 
regulation at 
Fansheiwe Dam and 
since 1965 at Wildwood 
Dam 
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TABLE 2. Comparison of Historical Means and Standard 

Deviations of Long-'Ierm and Short-Term Periods 





Long-Term 


Short-Term 


STATION 
02GD001 


HISTORICAL PERIOD 
(1516-13CG) 


riSTORICAL PERIOD 
(1923-1930) 




Mean 4 74 


516 




Std.Dev. 133 


88 


02GD003 


(1916-33,1936-44) 


(1923-1930) 




Mean 535 


586 




Std.Dev. 136 


127 



TABLE 3, 



Comparison of Monthly Means 



Means of Historical Record 



Station 


Jan, 


Feb. 


Mar, 


Apr . 


May 


Jtine 


July 


Aug* 


Sept. 


Oct. 


Nov. 


Dec. 




02GE002 


1333 


1558 


4116 


2561 


1162 


4 52 


386 


285 


325 


457 


1183 


1181 




02GD001 


608 


710 


1623 


1003 


461 


221 


171 


127 


146 


193 


Am 


472 




02GD003 


582 


730 


2157 


1188 


513 


176 


17 


101 


106 


169 


574 


572 












Means of 


Synthet 


ic Data 












I 

00 


02GE002 


988 


1323 


3926 


2586 


1317 


477 


363 


349 


273 


330 


98 3 


1154 


1. 


02GD001 


44 6 


630 


1556 


97 9 


490 


219 


164 


142 


131 


143 


397 


4 59 




02GD003 


423 


585 


2015 


1X87 


596 


193 


148 


135 


il 


95 


419 


531 





TABLE 4. 



Comparison of Standard Deviations of Monthly Ileans About Long-lerm Means 



Standard Deviations of Historical Record 



Station 


Jan. 


Feb. 


Mar - 


Apr. 


Ma^ 


June 


July 


Aug . 


Sept . 


Oct. 


Nov . 


Dec. 


02GE002 


113 9 


1694 


1261 


1494 


783 


272 


457 


199 


385 


534 


1403 


1081 


02GD001 


512 


697 


427 


503 


282 


143 


173 


69 


161 


233 


543 


408 


02GD003 


504 


872 


921 


735 


390 


118 


251 


105 


151 


251 


715 


541 



Standard Deviations of Synthetic Data 



a2GE002 


723 


1248 


1195 


2019 


1315 


339 


329 


297 


296 


262 


753 


1311 


02GD001 


294 


617 


442 


563 


352 


140 


117 


83 


125 


101 


276 


494 


02GD003 


415 


585 


889 


977 


644 


165 


179 


151 


102 


90 


341 


630 



TABLE 5 , Comparison of Skew Coefficients 

Skew Coefficients of Historical Record 

Station Jan. Feb. Mar. -Apr. May June July Aug. Sept, Oct, Nov. Dec. 

02GE002 0.7950 1.5860 0.6313 0.2592 0.0723 0.9998 1.7505 0.6322 1.9758 1.7421 1.3494 0.5168 

02GD001 0.8954 1.2035 -0.3580 0.3448 0.4154 1.5308 1.7797 0.4459 2.1597 2.0048 1.6457 0.5481 

02GD003 0.4463 1.8033 1.2220 0.1338 0.2237 0.7532 1.8548 1.3162 2.0903 1.6544 1.1365 0.4815 

I 

o 

Skew Coefficients of Synthetic Data 

02GE002 1.1181 1.4219 1.0512 2.0727 2.4423 1.5973 2.1035 1.6222 4.0475 2,1101 1.2827 2.0977 
02GD001 1.0666 1.4376 0.9672 1.4250 1.4523 1.2115 1.7305 1.0027 4.1328 1.5460 1.3887 2.2929 
02GD003 1.5069 1.6330 1.7258 1.8524 2.1349 2.1493 3.3216 1.9785 4.3208 1.9102 1.0009 2.1016 
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TABLE 6. Comparison of the Means, Standard Deviations and 
Coefficients of Skew for all Months 



Historical Record 



Station 
02GE002 
02GD001 
02GD003 
Total Skew 1.6490 



Mean 


Standard 
Deviations 

1437 


Coeff 
of 


■icients 
; Skew 


1249 


1. 


5571 


516 


558 


1. 


4245 


586 


760 


1. 


9655 



Synthetic Data 

02GE002 1172 1439 2.2721 

02GD001 479 533 1.8749 

02GD003 534 746 2.5572 

Total Skew 2.2347 
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TABLE 7 . Historical Correlation and Model Coefficieiit Matrices 

Mq - The Cross- Correlation Matrix (Variance-Covariance Matrix) 

02GE002 1.000 0.967 0.970 

02GD001 0.967 1.000 0.922 

02GD003 0.970 0.922 1.000 



M_, - The Lag-One Correlation Matrix (Lag-One Covciriance Matrix) 

02GE002 
02GD001 
02GD003 



A Matrix 

1.112 -0.221 '^ -0.523 

0.715 0.145 -0.473 

0.823 -0.267 _ -0.186 

B Matrix 

0.912 

0.877 0.237 

0.8 91 -0.058 0.223 



0.391 


0.373 


0.352 


0.396 


0.400 


0.353 


0.385 


0.358 


0.366 
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TABLE 8. Simulated Correlation Matrices 
50 Years of Synthetic Data 



M-. - The Cross-Correlation Matrix 



c:Ta^pJ^^] 


O^HFon? 


02GD001 


02GD003 


02GE002 


1.000 


0.972 


0.980 


02GD001 


0.972 


1.000 


0.923 


02GD003 


0.980 


0.923 


1.000 



M 



_, - The Lag-One Correlation Matrix 



02GE002 


0.419 


0.385 


02GD001 


0.426 


0.410 


02GD003 


0.403 


0.366 



0.384 
0.394 
0.385 
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TABLE 9. Comparison of Cross-Correlation Matrices on a Monthly Basis 



MONTH 



HISTORICAL DATA 



STATION 1 



January 



1* 
2 

3 



1.000 
0.998 
0,389 



0. 998 
1.000 
0.985 



0.989 
0.985 
1.000 



April 



1.000 
0.990 
0.988 



0.990 
1.000 
0.964 



988 
964 
000 



SYNTHETIC DATA 

ill 



1.000 
0.949 
0.912 



February 1 1.000 0,985 0.993 1.000 

2 0.985 1.000 0.957 0.980 

3 0.993 0.957 1.000 0.948 

March 1 1.000 0.904 0.916 1.000 

2 0.904 1.000 0.727 0.972 

3 0.916 0.727 1.000 0.961 



1.000 
0.959 
0.978 



0.949 
1.000 
0.840 

0.980 
1.000 
0.913 

0.972 
1.000 
0.919 

0.952 
1.000 
0,919 



0.912 
0.840 
1.000 

0.94 8 
0.913 
1.000 

0.961 
0.919 
1.000 

0. 978 
0.919 
1.000 



May 



1.000 
0.916 
0.984 



0.916 
1.000 
0.867 



984 
867 



1.000 



1.000 
0.952 
0.981 



0.952 
1.000 
0.932 



0.981 
0,932 
1.000 



June 



July 



1.000 
0.917 
0.841 

1.000 
0.999 
0.995 



917 
000 
584 



0.999 
1.000 
0.992 



0.841 
0.584 
1.000 

0.995 
0.992 
1.000 



1.000 
0.971 
0.970 

1.000 
0.959 
0.952 



0.971 
1.000 
0.908 

0.959 
1.000 
0.905 



0.9 70 
0. 908 
1,000 

0.952 
0.905 
1.000 



August 



1.000 
0.960 
0.921 



0.960 
1.000 
0.819 



0.921 
0.819 
1.000 



1.000 
0.958 
0.984 



0.958 
1.000 
0.922 



0.984 
0.922 
1.000 



September 1 
2 



October 



1.000 
0.993 
0.999 

1.000 
0.988 
0.994 



0.993 
1.000 
0.995 

0.988 
1.000 
0.9 78 



0.994 
0.995 
1.000 

0.994 
0.978 
1.000 



1.000 
0.986 
0.986 

1.900 
0.938 
0.942 



0.986 
1.000 
0.967 

0.938 
1.000 
0.852 



0.986 
0.967 
1,000 

0.942 
0.852 
1.000 



November 



1.000 
0.984 
0.995 



0.984 
1,000 
0.963 



995 
963 
000 



000 
896 
900 



0.896 
1.000 
0.755 



0.900 
0.755 
1.000 



December 



1.000 
0.985 
0.991 



0.985 
1.000 
0.953 



991 
953 
000 



000 
960 



0.981 



0.960 
1.000 
0,961 



0.981 
0.961 
1.000 



* NOTE: i:i the final generation of the data on the computer tapes 
Station 1 = 02GE002 

2 = 02GD001 

3 = 02GD003 
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T?vBLE 10 . ConTparison of Auto- Correlation Coefficient s 



HISTORICAL RECORD 



SYNTHETIC DATA 







Station 






Station 






ag 


1 


2 


1 


1 


1 




i 


§ 


1.00 3 


1.000 


1.000 


1.000 


l.OO'O 


1. 


000 


% 


0.406 


0.435 


0.338 


0.447 


0.477 


0. 


413 


1 


0.012 


O.OIS 


0,010 


0.101 


0.136 





.063 


f 


-0.113 


-0.12 3 


-0-077 


-0.03^ 


-0.040 


-0. 


056 


4 


-0.107 


-0.159 


-0.0 61 


-0.127 


-0,154 


-0 


.131 


1 


-0.2 30 


-0.2 70 


-0.221 


-0.216 


-0.277 


-0 


.215 


6 


-0.302 


-0.322 


-0.30 4 


-0.290 


-0.3 33 


-0. 


277 


f 


-0.241 


-0.272 


-0.2 33 


-0.24 8 


-0.290 


-0. 


.234 


i 


-0.144 


-0.172 


-0.123 


-0.139 


-0.170 


-0 . 


124 


i 


-0.061 


-0.056 


-0.066 


-0.085 


-0.083 


-0 


.076 


i 


0.082 


0.096 


0.071 


0.030 


0.055 


0. 


.022 



Station 1 = 02GE002 

2 = 02GD0 01 

3 = 02GD003 



TABLE 11. Results of Spectral Analysis 



Historical Data 
Station 

Cycle/Yr , 02GE002 02GD001 02GD003 



Synthetic Data 



0.0 


81 


0.2 


49 


0,4 


64 


0,6 


46 


0,8 


206 


1.0 


329 


1.2 


134 


1.4 


53 


1,6 


66 


1.8 


89 


2.0 


105 


2.2 


72 


2.4 


67 


2.6 


55 


2.8 


89 


3.0 


89 


3,2 


46 


3,4 


42 



4 
5 

9 
7 

3 3 
54 
22 

i 

11 
15 
15 
10 

ID 

i 

12 

12 

6 

i 



10 
14 
17 
12 
55 
87 
34 
12 
14 
20 
28 
19 
17 
14 
27 
29 
15 
13 







Station 




Cycle/Yr. 


02GE002 


02GD001 


02GD003 


0.0 


SO 


10 


19 


0.2 


m 


8 


15 


0.4 


71 


10 


18 


0.6 


$$. 


9 


16 


0.8 


232 


35 


58 


1.0 


349 


54 


90 


1.2 


120 


18 


32 


1.4 


51 


7 


15 


l.G 


Sfi 


8 


17 


1.8 


89 


12 


25 


2.0 


116 


15 


31 


2.2 


55 


6 


15 


2.4 


35 


4 


10 


2.6 


m 


5 


11 


2.8; 


ii 


11 


25 


3.0 


ii 


11 


25 


3.2 


II 


4 


i 


3,4 


24 


3 


7 
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TABLE 12 . xont for Long-lerm Persistence (Hurst':: K ) 

rtation Historic Record Synthetic Data 

02GE002 0.574 0.611 

02GD001 0.609 0.618 

02GD003 0.633 0.603 



TABLE 13. Comparison of the Number of Changes in Quartile 
per 50 Years 



Station Historic Record Synthetic Data 

02GE002 250 255 

02GD001 244 246 

02GD003 238 259 



TABLE 14, Connarison of the Number of Changes in the 
Direction of Movement 



Station Historic Record Synthetic Data 
02GE002 275 303 

02GD001 287 296 

02GD00 3 2 87 312 



Table 15. Comparison of Average Means and Standard Deviations of Daily Streamflows 



Historical Period 
1923-1930 



Historical Period 
1956-1970 



Station 



02GE002 
02GD001 
02GD003 

02GE002 
02GD001 
02GD003 

02GE002 
02GD001 
02GD003 

02 GEO 02 
02GD001 
2GD00 3 

02GE002 
02GD001 
02GD003 

02GE002 
02GD001 
02GD003 

02GE002 
02GD001 
02GD003 



Month 



May 



June 



July 



Aug. 



Sept, 



Oct. 



Nov. 



Mean 


Std.Dev 


1162 


1036 


461 


329 


514 


144 


452 


302 


222 


147 


176 


Ifti 


387 


382 


172 


169 


170 


304 


285 


205 


127 


75 


101 


131 


326 


241 


146 


114 


106 


118 


457 


213 


193 


90 


169 


154 


1183 


844 


460 


325 


574 


593 



Mean 


Std.Dev 


1162 


886 


481 


294 


548 


497 


499 


370 


224 


134 


199 


195 


331 


172 


170 


80 


10 8 


81 


366 


2 31 


195 


143 


119 


92 


346 


196 


173 


99 


117 


80 


479 


394 


183 


83 


216 


253 


1019 


80 8 


325 


202 


516 


537 



Synthetic 


Data 


50 yea 


irs 


Mean 


Std.Dev 


1322 


896 


489 


352 


596 


475 


478 


338 


218 


140 


193 


165 


367 


326 


164 


llf 


147 


178 


349 


217 


142 


83 


135 


121 


274 


246 


131 


125 


82 


102 


430 


262 


182 


101 


165 


m 


1000 


762 


397 


376 


519 


541 



I 

M 
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TABLE 16. ComparisonB of the Coefficients of Skew of Daily Streamf lows 



Historical Period 
1923-1930 



Station 


Month 


02 GEO 02 
02GD001 
02GD003 


May 


02GE002 
2GD0 01 
02GD003 


Jime 


02 GEO 02 
2GD001 
02GD003 


July 


02 GEO 02 
02GD001 
02GD00 3 


hug. 


02 GEO 02 
02GD001 
02GD003 


Sept. 


02GE002 
02GD001 
02GD003 


Oct. 


02 GEO 02 
02GD001 
02GD003 


Nov, 



Skew 

1570 
7846 
3 832 

2502 
1005 
1547 

4676 
3196 
9211 

9345 
0688 
5945 

90 58 
1078 
5542 



4. 
3 



351 
2362 
5018 



H i s to r ic al Pe r io d 
1956-1970 

Skew 

6.628 3 
5.7123 
6.52 80 

3.7458 
4.0574 
2925 



3.6120 
2 a 74 88 
4.2951 



4 
3 

5 

7 
9 
4 

7 

10 

2 

7 
3 
8 

2 
2 



0811 
3525 
825 8 

5640 
3466 
9909 

310 8 
2399 
9894 

3646 
72 89 
6790 

964 5 
2322 
3.4186 



Synthetic Data 
50 years 

Skew 

4.8186 
2.7101 
5.7723 

2.2222 
2.0886 
2.6604 

2 ,6278 
2.3131 
3.8547 

2.0156 
1,3406 
2.9982 

4.7033 
4.39 39 
5.1693 

2.8319 
2.1656 
2.7428 

1,9571 
2.1097 
2.2749 



TABLE 17. Coitiparison of Long-Term Period Means, May to November, Standard Deviations 
of Period Means About Long-Term Mean and Skews of Same 



Station 
02GE002 
02GD001 
02GD003 



Historical Period 
19 2 3-19 3Q 

Mean St, Dev. Skew 

607 1157 5.9745 

254 417 4.7010 

258 655 7,0134 



Historical Period 
1956-1970 

Mean St,Dev. Skew 

600 995 7.7853 

254 343 7.4096 

260 576 7.9902 



Synthetic Data 
50 years 



fean 


St.Oev. 


Skew 




589 


10 87 


6.8440 




240 


4 81 


3.6716 


1 
o 


239 


551 


8.4569 


1 
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Table 18. ?^nual Recorded Minimum Daily Discharges in CFS 



Year 

1916 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 



Minimum Daily Discharge 
Station: 02GE002 Station: 02GD001 



Station: 02GD00 3 



o 



d) 



145 July 9 
109 Nov. 11 
59.0 Sept. 5 
95.0 July 17 
130 Sept. 19 
162 Aug. 16 
106 Sept. 1 
60.0 Oct. 4 



199 
97, 
53. 
136 
140 
106 
79. 
93. 
60. 
56. 
102 
285 
43. 
107 
185 



Aug. 3 
June 20 
June 11 
Aug. 4 
Sept. 13 
Oct. 20 
July 14 
8 Aug . 2 6 
5 Oct. 30 
July 24 
July 12 
Sept . 6 
3 May 30 
Sept. 25 

Sept. 13 



en 



38.0 Nov. 6 




14.0 


Sept. 6 


38.0 Feb. 12 




15.0 


Sept. 23 


2 • Aug , 5 




20.0 


Aug. 7 


32.0 Aug. 11 




2.0 


Sept. 10 


4 8.0 Aug, 9 




26.0 


Sept. 2 7 


67.0 July 7 




24.0 


Aug. 2 


58.0 Sept. 5 




12.0 


Sept. 6 


37.0 Aug. 6 




16.0 


July 31 


58.0 Sept. 1 




19.0 


Sept. 16 


42.0 Sept. 3 




8.0 


Sept. 2 


40,0 July 26 




18.0 


July 20 


64.0 Sept. 26 




28.0 


Sept. 20 


80.0 Sept. 10 




34.0 


Sept. 10 


53.0 Sept. 9 




27.0 


Aug. 2 7 


50.0 July 21 




6.0 


Aug. 6 


57.0 June 24 




25.0 


Sept. 13 


79.0 Aug. 1 




35.0 


Aug. 2 


55.0 Nov. 6 




31.0 


Sept. 7 


50.0 July 2 






— 


31,0 Nov. 4 




^ 


-^. 


10.0 June 26 




7.0 


Aug. 18 


36.0 Aug. 1 




16.0 


Aug. 2 


35.0 Dec. 28 




17.0 


Dec. 2 


45.0 Oct. 7 




14.0 


July 25 


57.0 Aug. 12 




42.0 


Aug. 7 


29.0 July 2 




12.0 


Aug . 9 


46.0 Sept. 6 




18.0 


Aug . 2 2 


54.0 Oct. 11 




32.0 


Oct. 12 


54.0 Aug. 16 






— 


49.0 Jan. 1 






SSi<>S 


63.0 Sept. 30 






mam. 


114 Aug. 12 






^- 


92.0 Sept. 5 






-- 


92.0 June 27 






-- 


107 June 11 






-_ 


114 Aug. 19 









78.0 Oct. 27 






— 


58.0 Aug. 30 






-«, 


76.0 Aug. 29 




5.0 


Aug . 2 4 


63.0 Oct. 2 




1.0 


July 29 


9 3.0 Jan. 2 3 




38.0 


Oct. 26 


54.0 Aug. 23 




38.0 


Aug . 2 3 


63.0 Aug. 19 




4.0 


Sept. 16 


59.0 Aug. 16 




19,0 


Aug . 1 3 


73.0 Oct. 17 


^ 


3.6 


Dec. 9 


55.0 Jan. 24 





0.60 


Nov. 14 


44.5 July 16 


fe 


1.6 ] 


F'eb. 2 


45.6 Oct. 14 


-d 


0.50 


July 13 


54.6 Jan . 1 


4-* 


2.1 i 


Sept. 2 2 


54.6 Aug. 30 


cd 


0.80 


June 2 3 


60.0 July 24 


iH 


0.60 


Jan. 16 


112 Sept. 20 


& 


1.4 ] 


Dec. 5 


91.0 Sept. 1 


56.0 


July 31 


88.0 Aug. 11 




23.0 


July 16 



98.5 Sept. 2 



10.8 Nov. 19 
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# GAUGING STATIONS USED FOR STREAMFLOW 
GENERATION 



Figure I : Thames River Watershed 
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Figure 2: Hydrograph for Thames River Station: 02GE002, Year 1930. 
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Figure 3 : Hydrograph for Thames River Station: 02GD00I, Year 1930. 
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Figure 5 : Distributions of Historical and Synttietic Minimum Daily Flows, 
St 1 i on : 2 G D 1 , South Bra nch Th a m i s R i vtr . 
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Figurs 6: Distribution of HistOfical and Synthetic Minimum Daily Flows, 
Station: 2GD005, Nortti Branch Thames River. 
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